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Abstract 

Two types of particles, A and B with their corresponding antiparticles, are 
defined in a one dimensional cyclic lattice with an odd number of sites. In each 
step of time evolution, each particle acts as a source for the polarization field 
of the other type of particle with nonlocal action but with an effect decreasing 
with the distance: A ^ ■ ■■ BBBBB ■ ■ ■ ; B ^ ■ ■ ■ AAAAA It is shown 
that the combined distribution of these particles obeys the time evolution of 
a free particle as given by quantum mechanics. 

PACS 03.65.Bz, 03.65.Ca 
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Typeset using REVT^ 



The modeling of physical reality by means of fictitious particles that move and react in 
a substrate of different geometrical structures has been a fruitful strategy that has extended 
our analysis capabilities beyond the domain associated with differential equations [|I|. The 
particles involved in these models are classical in the sense that they are given precise lo- 
cation and velocity. This is clearly inadequate for the modeling of quantum systems that 
require, not only the indeterminacies imposed by Heisenberg's principle, but also nonlocal 
correlations between commuting observables, suggested by the Einstein Podosky Rosen ar- 
gument PI and empirically established in the violation of Bell inequalities 0,^. However 
this does not forbids the modeling of quantum systems if we do not identify the particles 
of the model with the quantum particles. It is possible, as will be seen in this work, to 
associate the real quantum particle to a combined distribution of two types of fictitious 
particles with nonlocal interaction. This simple example model can be trivially extended to 
higher dimensions of space and to a higher number of noninteracting quantum particles and 
it provides a new point of view to study the peculiarities of quantum mechanics. 

Let us assume a one dimensional lattice with sites on a circle and lattice constant a. 
We assume to be an odd integer. The reason for this restriction will become clear later. 
The inclusion of even values for would introduce unwanted complications in the model. 
Each site can be occupied by any number of particles of type A, B or by their corresponding 
antiparticles A, B. Particles and antiparticles of the same type annihilate in each site of 
the lattice leaving only the remaining excess of particles or antiparticles of both types A 
and B. At each time step, t — t + 1, corresponding to a time evolution by a small amount 
r, the particles of type A create antiparticles B in the same site, particles B in the first 
neighboring sites, B in the second neighboring sites and so on. In a similar way, particles 
B create particles A and A. 

A — B B B B B ■■■ 

B — y ■■■ A A A A A ■■■ . (1) 

The same reactions occur exchanging particles and antiparticles. This creation process 
extends to the right and left of each site up to the two opposing sites in the circle. Since 
A^ is odd, in these two sites particles of the same sign (either particles or antiparticles) 
are created. The number of particles or antiparticles created decreases with the distance d 
roughly like l/d"^ for a distribution of particles confined in a small region within a large lattice 
as will be precisely stated later. Before we write the master equation for the time evolution, 
we can notice some qualitative features of the process. It is easy to see that the process has 
diffusion. If we start, for instance, with some number of A particles in one site, after two 
time steps, some A antiparticles have been created at the same site reducing the number of 
A particles, but also some A particles appear in the first neighboring sites. The net effect is 
diffusion. It is less obvious that, even though the process has left-right symmetry, we may 
also have drift to the right or to the left. In order to see how this is possible we notice that 
A particles reject B particles from the site because B are created there, whereas B particles 
attract neighboring A particles to his site. Therefore if we have an asymmetric configuration 
like AB, the center of the combined distribution will move towards B. The drift direction 
and velocity is then encoded in the shape and relative distribution of both types of particles. 
We will see that, although the distribution of particles are widely distorted after few time 
steps, the drift direction and velocity remain invariant. 
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A convenient way to label the sites of the circular lattice is by an index s running from 
—L to L. Since the number of sites = 2L + 1 is odd, the index s will be integer. It is 
of course irrelevant which site has the label s = 0. Let as{t) and bs{t) be the number of 
particles of type A and B respectively at the site s at time t, normalized in a way that will 
be specified later (anyway the master equation is independent of the normalization). When 
as{t) or bs{t) take negative values they denote the number of antiparticles. At a particular 
site of the lattice, the number of particles change as particles or antiparticles are created in 
it by the particles in other sites. The time evolution of the process is then defined by the 
equations 



as{t + 1) = as{t) + E hs+d\{t)F{d) 

d=-L 
L 

hs{t + 1) = hs{t) - rg'' E , (2) 

d=-L 

where the square brackets in the index, \s + rf] , denotes "modulo A^" , that is, with a value 
in the closed interval [— g is related to the lattice constant ahy g = {2T[)/{Na) (it 
corresponds to the reciprocal lattice constant); r is a time scale small enough to make 
rg'^N'^ -C 1, and the function of the distance F{d) is defined as 

For later use we define a similar function G{d) as: 

G(d) = - E A;e*^'='^ = | 2siL7-V) if = ±1, ±2, ±2L _ 
^ ' N \ if d = ^ ' 

The alternating sign in the definition of F{d) corresponds to the fact that particles and 
antiparticles are created at alternating sites, and the different sign in Eq.2 is due to the 
difference in the role of particle and antiparticle in Relation 1. If the particles are confined 
in a small region within a large lattice, the main contribution in the sums of Eq.2 comes 
from terms with distance \d\ ^ A^. In this limit we have |-F((i)| ~ l/rf^ as mentioned before. 

The number of /I or 5 particles are not conserved in the time evolution. Neither is the 
sum of particles conserved. A quantity that is conserved in the time evolution of the process 
is the sum of the square of the number of A particles (or antiparticles) plus the sum of the 
square of the number of B particles (or antiparticles). This conserved quantity can be used 
for normalization and can be given a physical meaning like energy density or probability 
density. It is therefore relevant to define a combined distribution a^{t) + bl{t), associated to 
this conserved quantity. We will see that the drift velocity of this combined distribution, 
given by 

{V) = 4gJ2as{tMt)G{s-r) , (5) 

s,r 

is also conserved in the time evolution. Given a distribution {as{t),bs{t)}, we can change 
the drift velocity by an amount v, without changing the shape of the combined distribution 
by means of the local transformation 



3 



a'gif) = as(t) cos{vas/2) — bs(t) sin(fas/2) 

b'^{t) = asit) sm{vas/2) + cos(ms/2) . (6) 

These features have been checked in a computer simulation of the process. A circular 
lattice with = 801 sites (L = 400) and with lattice constant a = 1 was chosen. Several 
shapes of initial distributions were tried: gaussian, uniform and random, with several widths 
and drift velocities. The time dependence of M{t) = Ss('2s(^)+^s(^)) of the drift velocity 
given in Eq.5 was studied. Taking a time step r = 10~^, we found that these quantities 
remain constant after t = 1000 time steps, with a relative variation less than 10~^ for the 
gaussian case, 4 x 10~^ for the uniform distribution, and 0.04 for the random distribution. 
For larger time steps, r = 0.005, these quantities remain constant (less than 1% relative 
variation) for the gaussian and uniform case at t = 1000 but the random case begins to 
show significant departure from constancy. At r = 0.010 only in the gaussian case these 
quantities remain constant (less than 0.1%). The time evolution of the shape of the combined 
distribution strongly reminds the time evolution of quantum mechanical wave packets. For 
instance, a gaussian distribution for A and B particles, modified by Eq.6 in order to have 
drift, will evolve increasing the width and drifting but maintaining the gaussian shape. A 
uniform distribution will develop side lobes in the evolution. A remarkable feature is that 
the process smoothens out the random fluctuations of an initial distribution. 

The resemblance of the process with quantum mechanics is striking. We will indeed 
show that the process here defined corresponds to a quantum mechanical free particle in 
a lattice. Let us define then an N dimensional Hilbert space spanned by a basis {y:>s} 
s = —L, —L + 1, ■ ■ ■ , L corresponding to the eigenvectors of the position operator X. Then, 
Xips = as(fs- In this finite dimensional Hilbert space, we can not define the momentum 
operator P by means of the usual commutation relation. The alternative way to define P is 
to choose first an unbiased basis {4>k}, 

0. = 4^ E e^^^V. , (7) 

V A 

and with it, we define the momentum by the spectral decomposition 

L 

k=-L 

= l^T.9ke^'^'^'-'^^s{Vr,-) ■ (8) 

k,s,r 

The momentum eigenvalues and the relative phases to build the basis {4>k} have been chosen 
such that P is the generator of translations. That is, with this choice, the operator Ua = 
exp{—iaP) is such that Ua^Ps = V^s+i- The translation is cychc at the border, Ua^PL = V^-l- 
If we had taken A^ even, the right hand side of this equation should have a minus sign. This 
would comphcate the model of Relation 1 introducing a change of sign at some appropriated 
places. In order to have a simple lattice model for the quantum free particle we prefer to 
restrict ourselves to odd values of A^. 

The state of a free quantum particle, given by 
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^{t)= x: cs{t)^s , (9) 

s=-L 

will change according to the time evolution operator (we set h = 2m = 1) 

Ut = ex.p{-iPH) . (10) 

Let us consider the evolution of the coefficients of the expansion given in Eq.9, in one step 
of discretized time: to — rt and ti — T{t + 1) with a small time scale r and t positive integer. 
We have 

L 

Csit+1) = Crit){iPs,UriPr) ■ (11) 

r=-L 

For r small enough such that t||P^|| <C 1, that is r -C {a/nY, the time evolution operator 
can be linearized and we obtain 

L 

Cs{t + 1) = Cs{t) -ir Cr{t){Vs, P'^Vr) ■ (12) 

r=-L 

Using Eq.8 we calculate the matrix element 

(^„PV) = /^ E fcVt'=(-). (13) 

-'^ k=-L 

We have then 

L 

Cs{t + 1) = Cs{t) - irg^ Yl Cr{t)F{s - r) . (14) 

r=-L 

Reordering the terms in the sum and using the "modulo N" notation, we get 

L 

c,{t + 1) = c,{t) - irg'' Y cis+d\{t)F'{d) ■ (15) 

d=-L 

Finally if we explicitly write the coefficients with real and imaginary part, Cs{t) — as{t) + 
ibs{t), we get Eq.2 above. 

We can here check that M{t) — Y^g |cs(i)P is conserved. 

M{t + 1) = M{t) + r'g' ^ Cr{tK{t) ^ F{r - s)F{s - u) . (16) 

r,u s 

The term linear in r vanishes because the symmetric function F appears multiplied by an 
anti-symmetric factor. We see here that the "derivative" {M{t + 1) — M(t))/T vanishes like 
T in agreement with the numerical simulation of the process. We can write the functions F 
in their summation representations and, performing the sum over s, we get 

Y F{r - s)F{s -u) = ^ Y fcVi^*^^'-") . (17) 

s ^ k=-L 
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This sum can be evaluated as was done in Eq.3 and 4 but we don't need it. In the hmit 
N d = r — u, that is, when the particles are confined in a small region of a large lattice, 
we get 

M{t + l) = Mit)+T'^(^2Y,Cr{t)c:{t)^-0^ . (18) 

A similar result is obtained for the drift velocity, proportional to the expectation value of 
P. 

{P)t^-igY.<{t)cr{t)G{s-r) . (19) 

s,r 

In terms of as{t) and bs{t) this equation becomes the Eq.5 above. Here again, considering 
the time evolution {P)t+i, the term linear in r vanishes because it contains J2s[F{u — s)G{s — 
r) — G{u — s)F{s — r)] which is zero as can be calculated with the summation representation 
of the functions F and G. We obtain then 

(P),+i = (P), - tr'g' Y: c:{t)c,{t) E ' s)G{s - r)F{r - v) . (20) 

u,v s,r 

showing that the drift velocity is constant to order r, that is, the "derivative" vanishes 
with r in agreement with the numerical simulation of the process. Finally, applying a boost 
transformation exp{iXv/2) to the state of Eq.9, we prove equation (6). 

The one dimensional lattice model here presented provides a simple representation for 
the position and momentum of a free quantum mechanical particle. In this model we require 
that should be odd. Let us sec what happens in the case where N is even. In this case, 
the model evolves according to Eq.2 with the summations running from —N/2 to and 
with the same function F{d) defined in Eq.3. Notice that this function vanishes at the 
extreme values of d, that is F{±N/2) — 0. This model can be interesting in itself but it 
is no longer equivalent to the quantum mechanical system. The connection is lost in the 
step from Eq.l4 to Eq.l5. For the cases when the argument s — r of the function F in 
Eq.l4 take values exceeding N/2, we should introduce a minus sign if we want to change the 
argument to d as in Eq.l5 (in the case N odd, no sign change is needed). The reason for this 
change can be traced to the change in sign produced by the translation operator when the 
site labeled by ±L is crossed as mentioned after Eq.8. It would be possible to include even 
values for but at the cost of complicating the model. For this we would have to change the 
rules of Relation 1 exchanging particles and antiparticles when we cross the site with label 
±L. These complications are unwanted and we prefer to accept the fact that position and 
momentum of a quantum mechanical particle can be easily modeled only with a cyclic lattice 
with an odd number of sites. In the case of a quantum particle confined in a very small 
region (say, 10 sites) of a very large lattice (say, close to one million sites) it doesn't matter 
whether A^ is even or odd for all times until, due to drift or diffusion, the distribution 
reaches the sites with label close to ±N/2. However for small lattices and for extended 
distributions it does matter, and only in the odd N case the model of Relation 1 describes a 
quantum mechanical particle. This is a further indication of the essential nonlocal character 
of quantum mechanics. There is another case in quantum mechanics where an even or 
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odd number of states has important cualitative consequences. This is in finite dimensional 
reahzations of angular momentum. Whereas intrinsic angular momentum, the spin, of a 
particle can have an even or odd number of states, the orbital angular momentum, arising 
from position and momentum, can only have a realization with an odd number of states. 

The model presented can be extended from the free particle to the case of a position 
dependent potential. The general structure of the process shown in Relation 1 remains 
unchanged but the function F of Eq.2 will not be given by Eq.3 but will have to be calculated 
from an appropriate time evolution operator. The process can also be extended to two or 
three space dimensions but with larger computer requirements for the numerical simulations. 

Since the advent of quantum mechanics, there have been numerous attempts to develop a 
classical image for quantum behavior. For the reasons already mentioned at the beginning, 
the attempts in terms of particles are doomed. The model here presented, suggests the 
possibility of a classical image in terms of two fields A and B where each field acts as 
a source for the polarization of the other. As happens with the electromagnetic fields, the 
energy, or whatever conserved quantity, is given by the sum of the square of both fields. The 
consideration of these field may provide a new point of view for studying the peculiarities 
of quantum mechanics. 

I would like to thank H. Martin and A. Daleo for discussions and comments. This work 
has been done with partial support from "Consejo Nacional de Investigaciones Cientfficas y 
Tecnicas" (CONICET), Argentina. 
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